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The Brookhaven E821 experiment has recently mea- 
sured 1^ the anomalous magnetic moment of the muon 
with a precision of ±16 x 10"^", improving by a factor of 
6 on the previous measurement at CERN |^ . As the full 
set of data will be analyzed, the experimental error bars 
are expected to decrease further, by at least an additional 
factor of 3. The discrepancy between the experimental 
value of Ref. ||] and certain theoretical estimates can 
be as large as 2.6ct. Many possibilities to explain this 
difference through the introduction of new, beyond the 
standard model, degrees of freedom have therefore been 
considered. It is however important to keep in mind that 
the theoretical estimates of include several contribu- 
tions involving the nonperturbative hadronic sector of 
the standard model. In this Letter, we shall focus on 
one of these contributions, the so-called hadronic light- 
by-light scattering, which arises from the lowest order 
contribution, in the fine-structure constant a, to the ma- 
trix element (e stands for the electron charge) 

{^l-{p')\{^e)J,m^^-{p)) ^ (-^e)^i(p')fp(p^p)u(p) (1) 

of the conserved light quark electromagnetic current, 
jp — (2M7pU — d'-fpd — sjps)/3. This matrix ele- 
ment is given by a two-loop integral involving the con- 
nected fourth rank hadronic vacuum polarization ten- 
sor np^Ap(gi, 92, 93)- The corresponding contribution 
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This hadronic light-by-light correction to has been 
studied by several authors in the past ^-||. Its value, 
and even its sign, has suffered several changes, but the 
latest evaluations came to values that agreed within the 
quoted theoretical error bars and that were negative. In 



particular, there exists a well defined contribution to 
^LbyL; had ^ dcuotcd by a}"p^'^'''^ , that arises upon restrict- 
ing flp^Apl?!, (72, 93) to its reducible one-pion exchange 
component and which represents about 70% of the total 
value of a^^^'^'- Very recently, two of us, reconsid- 
ering this evaluation of a^i^^^''^^ , found a positive result 
0, but which, in absolute value, agreed with the pre- 
vious numerical values whenever comparison was possi- 
ble. Resolving this sign ambiguity certainly represents 
a major issue, since the result of |^ reduces the the- 
ory vs experiment discrepancy to less than 1.5a. This 
Letter aims at providing an argument in favor of a pos- 
itive sign for a^^^'^ which, to a large extent, does not 



rely on the methodology followed in Ref. 



The off- 



shell pion-photon-photon form factor T^o^^,^* constitutes 
an important ingredient for the evaluation of a}"^^^''^" . 
For a constant form factor, fixed by the Wess-Zumino- 
Witten (WZW) term of QCD 0, which reproduces the 
Adler ||ll|, BcU-Jackiw jl^ anomaly, the corresponding 
two-loop integral for a^^^^'-'^ diverges like {a/i:YC Iv? A, 
where A stands for an ultraviolet cutoff jl3j. We shall 
be interested in the determination of the coefficient C of 
this log-squared divergence, so that we may then compare 
it to the value that can be extracted from the analysis 
of Ref. 1^. In order to achieve this goal, we shall use 
a renormalization group argument, along the lines dis- 
cussed in Ref. |jlj] , within the framework of the effective 
low-energy field theory of the standard model. Before 
doing so, we first describe the relevant features of this 
effective theory. 

The low-energy degrees of freedom of the standard 
model involve the pseudoscalar mesons, the light lep- 
tons, and the photon. The interactions between these 
degrees of freedom are constrained by the symmetries of 
the standard model, like chiral symmetry or C/(l) gauge 
invariance. This effective theory will provide a good 
description of physical processes as long as the energy 
scales involved are much lower than a typical hadronic 
scale Kfi ^ 1 GeV. Since the leptons obey a first order 
equation of motion, the chiral counting has to be suit- 
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ably adapted. The situation differs from the case of chi- 
ral perturbation theory in the presence of baryons due 
to the fact that the mass m of the light leptons does 
not introduce a scale of the order of A^- Rather, the 
counting m ^ C(p) arises quite naturally. Furthermore, 
we shall also count fermion bilinears as quantities of or- 
der p, as is the case in the sector of the effective the- 
ory which describes the semileptonic decays of the pseu- 
doscalar mesons (see also Ref. for a discussion 
of this aspect in a different context). Finally, the elec- 
tric charge is also counted as a quantity of order p [ p7[ . 
The lowest order term in this extended chiral expansion 
then starts at order p^, with (for the notation, see Refs. 
RPIPI): 



+ ^{{d^U+df,U)+2Bo{M+U + U+M)). (3) 

This lowest order Lagrangian involves the Maxwell La- 
grangian for the photon, the tree-level minimal coupling 
of leptons and photons, as well as the lowest order chi- 
ral Lagrangian for the mesons, together with their cou- 
plings to the photons. It gives the (finite) lowest order, 
0{p^), contribution to a^^^' obtained upon restricting 
np,yAp(<Zi, 92, 93) to its lowest order approximation, con- 
sisting of a charged pion loop with pointlike electromag- 
netic vertices. At 0{p'^) in the Lagrangian, we find the 
usual counterterms involving the low-energy constants Li 
@, Ki and X, These counterterms absorb 

the divergences due to loops with virtual pions, photons, 
and leptons. For instance, Lg is needed in order to make 
the slope of the pion form factor, generated by a charged 
pion loop, finite. Note that the pion loop contribution to 
^fj,i^Xp{qi, 92, 93) with these nonpointlike vertices, gives a 
next-to-leading, 0(p®), correction to aj;''^'"' which is no 
longer finite, and thus requires a corresponding countert- 
erm from £^^^^ . In addition, we also encounter in jC'-^^ the 
counterterms which result from the divergent loop graphs 
involving virtual leptons and photons, such as the 0{a) 
renormalizations of e, m and of the lepton wave function. 
These will play no role in the present discussion. Finally, 
the WZW Lagrangian of the odd intrinsic parity sector 
also occurs in C'^^\ For our present purpose, we need 
only to retain a piece of the latter (our conventions are 
as follows: 75 = ij^^f-^j'^^^, £0123 = +1) 
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This term is responsible for the lowest order, and diver- 
gent, contribution to the pion-pole correction a^^^''^ , 



which thus starts at order p* and corresponds to the dia- 
grams (a), (b), and (c) of Fig. 1. At order p^ in the effec- 
tive field theory, we encounter, for instance, divergent 
loops involving the particular WZW vertex from 

£(4) 

shown in Eq. (^ and a virtual fermion line, like precisely 
the triangular subgraphs in the two first graphs of Fig. 1. 
These divergences have to be canceled by an appropriate 
set of counterterms which, according to Ref. pl| , reads 



£(6) 
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Here also, there may be other terms, indicated by the 
ellipsis, but which are of no relevance for the discussion 
concerning a^^^'^ . Under renormalization, the terms 
we have shown absorb the divergence of the triangular 
subgraphs in the diagrams (a) and (b) of Fig. 1. They 
therefore also contribute to a^^^''^° at order p^, through 
the diagrams (d) and (e) of Fig. 1 . We shall come back to 
this important issue below. Let us first conclude this brief 
description of the effective theory by mentioning that 
the graphs (a) and (b) also have an overall divergence, 
which is absorbed by an appropriate set of counterterms 
in C^^^\ the effective Lagrangian at the ©(p^") level. 
The corresponding tree-level contribution to a^^^'"^ is 
represented by the diagram (f) in Fig. 1. The detailed 
structure of C^-^*^^ will, however, not be needed. 
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FIG. 1. The graphs contributing to " ' at lowest order 
in the effective field theory. The two-loop graphs (a), (b), and 

(c) contain two WZW vertices from The one- loop graphs 

(d) and (e) have an insertion of one vertex (■) of C^^\ Eq. 
(^). Finally, the tree-level graph (f) stems from the ©(p^^) 
counterterms. Note that the diagram (c) is actually finite. 

We next return to the original purpose of this Letter, 
the determination of the two-loop leading logarithm aris- 
ing from the sum of the diagrams (a) and (b) of Fig. 1. 
Since the diagram (c) is finite, it is of no relevance for 
the present discussion. In the context of a renormaliz- 
able field theory, such a question would be most naturally 
addressed by means of the renormalization group equa- 
tions. As stressed by Weinberg the renormalization 
group remains a useful concept even in the case of a non- 
renormalizable local effective field theory. Within the 
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effective theory framework, the expression of the renor- 
mahzed contribution to arising from the graphs of 
Fig. 1 takes the general form 



LbyL;7r 



H 
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Here denotes the arbitrary subtraction scale introduced 
by the renormalization procedure. The contributions of 
the two- loop graphs (a), (b), and (c) of Fig. 1 are given by 
the first term in this expression. The next term describes 
the contributions from the one- loop graphs (d) and (e), 
with an insertion of the renormalized 0(p®) counterterm 
X = (— l/4)(xi -I- X2), while the last term, collec- 
tively stands for the renormalized tree-level contributions 
from the ©(p^") effective Lagrangian. The dependence 
on the subtraction scale fj, in the two- and one-loop func- 
tions reads 



H 



= J2 ^p^^" 
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The (dimensionless) coefficients hp and jq do not depend 
on /i, but are functions of the ratios M^o/m and FTj/m. 
They therefore satisfy Vhp = 0, Vjq = 0, where 
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Being a physical quantity, a^^^^''^ obeys the condition 
^ji{da^^''^'-'^° / d^i) = 0, which gives 



In 
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= 0, 
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where we have introduced 7^^ = ^[dx{^Jt) / d^\. Acting 
with the operator T) on this equation and using the fact 
that Vxi^j) = 0,T>k{^) = 0, one finds that the term 
proportional to ln(TO//i) has to vanish separately, i.e., 
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Therefore, in order to compute h2 = {a/iry^C, we have 
to extract the dependence on h\{m/ ^l) from the one- loop 
graphs of Fig. 1. This is our next step. 
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This integral diverges logarithmically. Introducing a cut- 
off A, the corresponding contribution to is obtained 
upon using (||) . The Dirac trace can easily be performed, 
leading to the result a^^^-'^ l(d)+(c) = ji ln(TO/A)-|-finite, 
with 
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It is interesting to note that the result for ji can also 
be obtained from the calculation of the coefficient of the 
ln(M^/TO'^) term in the two-loop electroweak contribu- 
tion to a^; see Ref. |22 . Next, we recall that the scale 
dependence of the constant x(/i) is already known pl| , p3t . 
Indeed, the same combination of xi ^^^d X2 arises in the 
0{a^) contribution to the tt" e+e^ amplitude, to- 
gether with the divergent one-loop graph, generated by 
the WZW vertex, and involving two virtual photons and 
one fermion line [the on-shell restriction of the triangu- 
lar subgraphs in diagrams (a) and (b)]. The result reads 
7^ = Nc, and thus with Eqs. (^ and ( p^ ) we obtain 



where 
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For F^ = 92.4 MeV, to = 105.66 MeV and Nc = 3, 
this gives C — 0.0248. We shall now discuss some conse- 
quences of this result as far as the sign of a}"^^^'-'^ obtained 
in Ref. ||^ is concerned. 

Our analysis thus tells us that the sum of the graphs 
(a) and (b) of Fig. 1 not only diverges, but behaves like 
{a/'KY'C \t? a as the ultraviolet cutoff A is sent to infinity. 
Actual calculations of a^^^^''^ [which, in the framework 
of the effective theory, amount to estimates of x(m) and of 
K(/i)] replace the pointlike WZW vertex by a form factor 
T-^a^-i-, which regularizes the ultraviolet behavior. The 
models of J-T^a^.^. that are introduced for this purpose in 
one way or another involve a hadronic scale Ajj, typically 
the p- meson mass My (or Ajj ~ 2toq, with a constituent 
quark mass toq.). As My 00, the form factor becomes 
constant, and one therefore has to recover the same be- 



The sum of the two graphs (d) and (e) of Fig. 1 reads havior, a]^"^'-'^ ~ {a/irfClT? My ■ To the best of our 
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knowledge, this aspect of a^^^^''^ has not been discussed 
in detail before. In particular, we know of no previous 
determination of the constant C. There is a reference to 
the In^ A behavior in [T^ , which, in our notation, states 
that a'i^''^-''°\wzw ~ {a/'Kf{m/'iTTF^{\T?K)/2, but it 
is not clear from the context of Ref. whether the sign 
was meant to be part of the estimate or not. 

In Ref. a^^^^'^ was expressed, for a certain class 
of form factors, in terms of a two-dimensional integral 
representation of the form (we omit here the small con- 
tribution from the finite diagram (c) in Fig. 1, which 
plays no role in the present discussion) 



The implications of these observations for the comparison 
between theory and experiment as far as the anomalous 
magnetic moment of the muon is concerned have already 
been discussed in Ref. ||]. In particular, the difference 
between the present experimental value of 0] and the 
standard model is reduced from 2.6(t to 1.5(7. The exten- 
sion of the present effective field theory analysis to the 
other contributions from hadronic light-by-light scatter- 
ing to the muon g — 2 will be discussed elsewhere. 

This work has been supported in part by Schweiz- 
erischer Nationalfonds and by TMR, EC Contract No. 
ERBFMRX-CT980169 (EURODAPHNE). 



dQi dQ2[f{Mv,QuQ2)wfiQi,Q2) 
Jo 

+g{Mv,Qi,Q2)wg{Mv,Qi,Q2)]- (16) 



The functions f{Mv,Qi,Q2) and g{Mv ,Qi,Q2) de- 
pend quadratically on J-^^o^,^,. The weight functions 
Wf{Qi,Q2) and Wg{Mv ,Qi,Q2) are both positive and 
peaked around Qi ^ Q2 ^ 0.5 GeV. In the case 
of the vector meson dominance (VMD) form factor, 
/^*^-° (My, Qi, (52) vanishes, and only the term involv- 
ing Wg{MY,Qi,Q2) contributes. The corresponding 
value obtained in ||] for a^^'-'^°\vMD upon numerical 
evaluation of the integral dlfl) differs only by its overall 
sign from the results obtained by previous authors [^,^ . 
Now, in Ref. it has been shown that (up to ~ InAfy 
or constant terms) 

poo poo 

lim / dQi dQ2g''^'''{Mv,Qi,Q2)wg{Mv,Qi,Q2) 

= Cln^My, (17) 

with a constant C that numerically agrees with the re- 
sult of Eq. (^5[), including its sign. A global sign error 
in the second contribution to (16) is thus excluded. On 
the other hand, for a constant form factor, normalized by 
the WZW term, f^^"^ {My , Qi, Q2) = [Nc / [Un^ F^)]^ 
and g^^^ {Mv,Qi,Q2) vanishes. In that case, the in- 
tegral involving Wf{Qi,Q2) diverges. In Ref. it was 
shown that this divergence is of the form (up to ~ In A 
or constant terms) 



lim / dQi dQ2r^"'{Mv,Qi,Q2)wfiQi,Q2) 

^Cln^A, (18) 

where, again, the constant C was found j^] to be positive 
and in agreement with Eq. dl^). We thus conclude that 
the overall sign of the first term in Eq. (|l^) also has to 
be correct. 

Therefore, we expect that a'^^^^''^ is a positive quan- 
tity, at least for the VMD type of parametrizations of 
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